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ABSTRACT 

It is widely thought that core-collapse supernovae (CCSNe), the explosions of massive stars following the 
collapse of the stars’ iron cores, is obtained due to energy deposition by neutrinos. So far, this scenario was not 
demonstrated from first principles. Kushnir and Katz (2014) have recently shown, by using one-dimensional 
simulations, that if the neutrinos failed to explode the star, a thermonuclear explosion of the outer shells is 
possible for some (tuned) initial profiles. However, the energy released was small and negligible amounts of 
ejected o6 Ni were obtained, implying that these one-dimensional collapse induced thermonuclear explosions 
(CITE) are unlikely to represent typical CCSNe. Here I provide evidence supporting a scenario in which the 
majority of CCSNe are the result of CITE. I use two-dimensional simulations to show that collapse of stars that 
include slowly (few percent of breakup) rotating ~ 0.1 — 10 M 0 shells of mixed helium-oxygen, leads to an 
ignition of a thermonuclear detonation wave that unbinds the stars’ outer layers. Simulations of massive stars 
with different properties show that CITE is a robust process, and results in explosions with kinetic energies 
in the range of 10 49 — 10 52 erg, and 56 Ni yields of up to ~ Mg, which are correlated, in agreement with 
observations for the majority of CCSNe. Stronger explosions are predicted from higher mass progenitors that 
leave more massive remnants, in contrast to the neutrino mechanism. Neutron stars are produced in weak 
(<; 10 51 erg) explosions, while strong (> 10 51 erg) explosions leave black hole remnants. 

Subject headings: hydrodynamics — methods: numerical — supernovae: general 


1. INTRODUCTION 

There is strong evidence that CCSNe are explosions of 
massive stars, in volvin g the collaps e of the stars’ iron cores 
dBurbidge et all [19571 Hirata et al.l 119871 ISmartd l2009h and 
ejection of the outer layers. It is widely thought that the 
~ 10 51 erg observed kinetic energy (KE) is obtained due to 
the deposition of a small fraction (~ 1%) of the gravita tional 
ener gy (~ 10 5 3 erg) released in neutrinos (see lBethall990t 
iBurrowsI 120131 . for reviews). So far, this scenario was not 
demonstrated from first principles. 

iBurbidge et all (119571) suggested a different mechanism for 
the explosion during core-collapse that does not involve the 
emitted neutrinos. In the proposed scenario, increased burn¬ 
ing rates due to adiabatic heating of the outer shells as 
they collapse lead to a thermonuclear expl osion (see also 
iHovle & Fowlerfl960l:lFowler & Hovlell 19641) . This collapse- 
induced thermonuclear explosion (CITE) naturally produces 
~ 10 51 erg from thermonuclear burning of ~ M 0 (gain of 
~ MeV/m „). A few one-dimensional studies suggested that 
CITE f ails ((Colgate & WhitefI%6tlW6oslev & Weaveill98l 

iBodenheimer & Wooslevlll983l) . It was recently shown, by 
using one-dimensional simulations, that CITE is actually pos¬ 
sible f or some (tu ned) initial density and composition profiles 
dKushnir & Katzl2014t) . which include shells of mixed helium 
and oxygen. An ignition of a thermonuclear detonation is ob¬ 
tained, ~ 10 seconds after the core collapsed, that unbinds 
the outer layers of the star. However, the energy released 
was small, < 10 5 ° erg, and negligible amounts of ejected 56 Ni 
were obtained, implying that these one-dimensional CITE are 
not typical CCSNe. 

Here I provide evidence supporting a scenario in which the 
majority of CCSNe are CITE of rotating massive stars. By 
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using two-dimensional simulations with a fully resolved igni¬ 
tion process, I show that collapse of stars that include slowly 
(few percent of breakup) rotating ~ 0.1 — 10 M 0 shells of He- 
O with densities of few x 10 3 g cm -3 , leads to an ignition of a 
thermonuclear detonation that unbinds the stars’ outer layers. 

Previous studies of the rmo nuclear expl osio ns of col¬ 
lapsing rotating stars (IBodenheimer & Wooslevi 119831 : 
iMacFadven & Wooslevl 119991) . without He-O mixtures, did 
not result in an ignition of a detonation. I show that the 
presence of He-O mixtures is a necessary condition for 
explosion, making the previous studies consistent with the 
present work. 

Section[2]describes the simulations. Section[3]demonstrates 
the consistency of CITE with the majority of CCSNe, includ¬ 
ing new predictions. 

2. NUMERICAL SIMULATIONS 

Section IT71 describes the initial pro files. Section [272] de¬ 
scribes the numerical tool. Section [2731 presents a simulation 
of a successful ther monu clear explosion that leads to a typical 
CCSN, and Section [2741 demonstrates that CITE is robust, by 
examining the sensitivity of the results to the assumed initial 
profile. 

2.1. Initial profiles 

Pre -collapse profiles are uncertain (see, e.g lSmith & Arnettl 
120141) . but nevertheless stellar evolution calculations of ro¬ 
tating m assive stars generally predict the existence of a He- 
O she ll dHirschi et al]l200i I200I lHirschj|l2007l lYusofetalJ 
120131) . Instead of relying on some specific calculated pro¬ 
files, I systematically examine the explosion’s properties as a 
function of the He-O shell properties. Profiles are composed 
as f ollows (this naram etrization is somewhat different from 
iKushnir & Katzl (20141) ): 

• The profile is composed of shells with constant entropy 
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(per unit mass), s, constant composition, and in hydro¬ 
static equilibrium. 

• 1.2 Mg is placed within r < 2 • 10 8 cm, representing a 
degenerate iron core. 

• The entropy of the He-O shell is s H e-o and the mass 

fractions are A Hc = 1 — Xq. At the base of the shell the 
enclosed mass is Mt> asc , the density is phase, and the ratio 
between the local burning time and the free-fall time, 
tb/tff, is where tb = e/Q , £ is the internal 

energy, Q is the thermonuclear energy production rate, 
and tff = r 3 / 2 / s/2GM{r). The inner radius of the 
shell, rbase, is set by these parameters, and the outer 
radius is placed at p = 10 3 g cm~ 3 . 

• Pure oxygen (helium) is placed below (above) the He-O 
shell, where the entropies are determined by requiring 
enclosed mass of 1.2 Mq within r < 2 • 10 8 cm and 
zero pressure or low temperature (10 5 K) at the profile’s 
fixed outer radius of 3 • 10 10 cm. 

Instead of specifying .s Hc -o, it is more convenient to specify 
M S heii, the mass of the He-O shell. For given Mb ase , Phase, and 
tbfi/tfffi, there is a maximal possible M s h e ii (obtained for a 
minimal possible sne-o), denoted b y Mshdi.max. _ 


Following previous studies ( 

Bodenheimer & Wooslevi 

fl983l iMacFadven & Wooslel 1999 

), the angular momentum 


is initially distributed such that / rot , the ratio of the centrifugal 
force to the component of gravitational force perpendicular 
to the rotation axis, is some constant / m t,o (0.05 for most 
simulations) throughout the profile, except for the following: 

* ,/rot — 0 at small radii, r ryot,o = -rmLOr*base ■ -riot.o — 
0.8 for most simulations. 

• /rot = 0 at large radii, and increases linearly with de¬ 
creasing radius between r = 2 • 10 10 cm and r = 
1 ■ 10 10 cm to /rot.o- This is done for numerical stability, 
and has a small effect on the results. 

The initial profiles are chosen to be in hydrostatic equilibrium 
without rotation. Rotation is added such that initiall y the pro¬ 
file is slightly out of equilibrium (see Section [2.4.2l >. 

2.2. Collapse Simulations 

The problem considered is axisymmetric, allowing the use 
of two-dimensional numerical simulations with high resolu¬ 
tion. I employ the FLASH4.0 code with thermonuclear burn¬ 
ing (Eulerian, adaptive mesh refinement; iDubev et al.1120091) 
using cylindrical coordinates (II, z) to calculate one quadrant, 
with a resolution (i.e. the minimal allowed cell size within the 
most resolved regions) of ~ few x 10 km. I use th e 13 isotop e 
alpha-chain reaction network of VULCAN 2D dLivnell 19931) . 
and a burning-limiter dKushnir et al.ll2013l) . 

Angular momentum was implemented similarly to previ - 
ous studies dLindner et al.l2010tlFernandez & Metzgerl2013l) . 
Specific angular momentum is included as a mass scalar quan¬ 
tity. The centrifugal force is included as part of the default 
treatment of “fictitious” forces that arise in curvilinear coor¬ 
dinates. The specific KE associated with rotation is added to 
the specific KE used in the regular scheme of FLASH4.0, as 
necessary when viscous forces are neglected. 

I assume that neutrinos escape with negligible effect on the 
outer layers (gravitational mass loss is only considered for the 


post-processing estimate of the remnant mass). Layers below 
the inner boundary, ri nner (5 • 10' cm for most simulations), are 
assumed to have already collapsed, and the pressure within 
this radius is held at zero throughout the simulation. Mass 
of material that (freely) flows through the boundary is added 
to the original collapsed mass and contributes to the gravita¬ 
tional field. 

At locations in the progenitor where T > 2 • 10 9 K, oxygen 
is replaced with silicon to prevent fast initial burning, and the 
region between r; nner and r = 2-10 8 cm is filled with s = Iks 
iron (satisfying hydrostatic equilibrium), which is prevented 
from burning. These regions collapse rapidly thorough the 
inner boundary and have negligible effect on the results. 

For simulations in which an explosion was obtained, the 
original resolution was reduced successively by factors of 2 
for each doubling of the shock radius, starting with the first 
resolution increase made at a chosen threshold radius r res (5 • 
10 9 cm for most simulations). 

2.3. CITE is possible - example of a successful explosion 
which leads to a typical CCSN 

A pre-collapseprofile which leads to a typical CCSN is 
shown in Figure [I] The dynamical evolution of the collapse, 
calculated with a resolution of ps 30 km, is shown in Figure^ 
Collapsing material that includes angular momentum reaches 
a centrifugal barrier at a radius smaller by ~ / rotj o of its initial 
radius, leading to an increased pressure and to the formation 
of a rotation-induced accretion shock (RIAS), seen in panel 
(a). The in-fall velocity of the shock-heated oxygen is sig¬ 
nificantly reduced, allowing synthesis of 56 Ni (black contour) 
within the hot and dense downstream conditions. The energy 
released in this thermonuclear processing (~ 3 • 10 5 ° erg un¬ 
til t = 28 s, panel (d)) is sufficient for an expansion of the 
synthesized material. The base of the He-O shell, seen in 
panel (a), is deformed during collapse from its initial spherical 
shape. The RIAS hits the He-O shell at t ~ 29 s which causes 
an ignition of a detonation (panel (b)). Because the shapes of 
the RIAS and the base of the He-O shell are different, the igni¬ 
tion position propagates from the symmetry axis towards the 
equator. The detonation propagates outward (panel (c)), pro¬ 
ducing thermonuclear energy at fewx 10 5 ° erg s -1 (panel (d)). 
The pressure built from accumulating thermonuclear energy, 
aided by the RIAS that acts as a piston, manages to halt the 
inward collapse and causes an expansion that leads to an out¬ 
ward motion. Once the detonation reaches outer layers with 
densities p < few x 10 3 g cm -3 it decays and transitions to a 
hydrodynamic shock which continues to propagate outwards. 

The resulting ejecta has a mass of ps 4.5 Mq (leaving 
pa 6.5 Mq of material with negative spherical radial veloc¬ 
ity, hereafter in-falling mass), KE sa 1.3 ■ 10 51 erg, and a 56 Ni 
mass, M n ;, of ~ 0.08 Mq. The obtained KE and M N ; are 
typically observed in CCSNe (see Section [3|. The KE of the 
ejecta may change slightly if a hydrogen envelope is added. 

2.3.1. Numerical convergence 

The KE and Mn; as a function of resolution are shown in 
panel (a) of Figure[3] The results are converged to ps 10% for 
resolution of « 30 km, used throughout the paper. 

The 56 Ni yield is estimated as the total mass of 56 Ni with 
positive (spherical) radial velocity, Mm(v r > 0), at the latest 
time with the original resolution. At this time Mm(v r > 0) 
is roughly constant with time, and a decrease in Mm (v r > 
0) was obtained as the resolution was reduced. In order to 
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verify that the decrease in Mm(v r > 0) is numerical, another 
calculation with r res = 10 10 cm was performed (x-symbols 
in panel (a) of Figure in which similar Mm was obtained. 
The reduction of the resolution inhibits accurate calculation 
of the asymptotic 56 Ni distribution in the ejecta. 

The KE and Mm as a function of rj nnel are shown in panel 
(b) of Figure [3] The smaller r; nner , the earlier the RIAS is 
launched. As can be seen the KE and Mm are not sensitive 
to the exact value of rj nner with the Mm slowly increasing for 
smaller ri nner . It would be desirable to decrease ri nner below 
the minimal value tried here (100 km), however the in-fall ve¬ 
locities at this radius are already « 0.4c, leading to tens of 
percent error for a non-relativistic calculation. In what fol¬ 
lows I employ rj nnel - = 500 km, keeping in mind that Mm is 
probably under-predicted by a factor of a few. 

2.4. CITE is robust - the full set of simulations 

In Sections [2.4. ll[2~4~4l I examine the sensitivity of the re¬ 
sults of Section |23lto z rnt n, / rot ,o> -^He, M she h, and 4,o/i//,o, 
while in Section 12.4.51 M base and p base are varied. The relevant 
figure for each section is indicated in the title. A list of the 
simulations is given in Tables [lEl 

2.4.1. Sensitivity to the inner boundary of the rotation zone (x w t, o) 
- panel (c) of Figure\3\ 

For 0.6 < a; ro t,o < 1 ignition occurs immediately after the 
interaction with the RIAS, similarly to panel (b) of Figure 0 
The smaller x rot .o, the earlier the RIAS is launched, leading 
to smaller KE (since the He-O shell is less compressed) but to 
higher Mm (because of the higher density of the oxygen). 

For :/: r() | 0 = 1 there is no RIAS inner to the He-O shell, 
and it ignites because of adiabatic compression during free- 
fall Jselfzigtirtirm),jiimilat:\y to one-dimensional simulations 
dKushnir & Katzll20l4[) . In fact, in a one-dimensional simula¬ 
tion of this profile, while the explosion fails, « 2.5 • 10 51 erg 
of thermonuclear energy is released, which is comparable to 
the * roti o = 1 case (« 2.9 • 10 51 erg). The one-dimensional 
case fails because only a small fraction of the thermonuclear 
energy is converted to outward motion and cannot overcome 
the binding energy of the star ([Kushnir & Katzll2014l) . A’ b i n ~ 
—9 • 10 5O erg (exterior to the base of the He-O shell, cor¬ 
rected for thermal energy). For :£ rol .o = 1, the RIAS formes 
inside the He-O shell after ignition, and reaches the detona¬ 
tion front while acting as a piston that increases the fraction 
of thermonuclear energy that is converted to outward motion 
(« 0.75 in this case, < 0.35 in the equivalent one-dimensional 
case). Small amounts of 56 Ni are synthesized, since the RIAS 
forms late, where only a small amount of high density ma¬ 
terial is present. As .x mL () is increased, a smaller fraction of 
the He-O shell is RIAS supported, and for x' roL o = 1.2 the 
explosion fails, similarly to the one-dimensional case. 

For xvot.o <0.7 more complicated dynamics are obtained, 
since the RIAS hits the He-O shell at early times, in which 
the density is too low for immediate ignition. However, later 
on ignition happens at the pole, where the induction time is 
shortest, and then the detonation slides towards the equator. 
The KE is smaller by a factor of a few compared to shock 
ignition, but the Mm is somewhat higher because of the high 
density material that is being shocked. 

2.4.2. Sensitivity to the rotation speed (f m t,o) - panel (d) of 
Figure\3\ 


The KE depends weakly on / ro t,o> as long as f mt o > 0.02. 
For smaller values, KE drops sharply, and no explosion is ob¬ 
tained for / ro t,o = 0.01, as no RIAS was launched at rele¬ 
vant times. The Mm decreases with increasing / rot .o, since 
the RIAS is launched at larger radii, leading to collapsing ma¬ 
terial with smaller density. This demonstrates that rotation is 
required for CITE. 

To test the effect of the artificial departure from hydrostatic 
equilibrium due to rotation, I preformed a set of simulations 
where the initial profiles were first relaxed to equilibrium be¬ 
fore the collapse. The relaxation was run for ~ 1000 s without 
inflow through ri nner and without burning. The results of such 
simulations (/ ro t,o = 0.05 and / roti0 = 0.1) change by <; 50% 
(x-symbols). 

2.4.3. Xne - panel (e) of Figure\3\ 

For low (high) values of A He the KE decreases since the 
energy content of the He-O mixture is decreasing (since there 
are not enough target nuclei for alpha capture and the energy 
content of the He-O mixture cannot be extracted efficiently), 
until the explosions fails for Ane = 0.1(0.9). This demon¬ 
strates the He-O mixture requirement for CITE. 56 Ni is pro¬ 
duced below the He-O shell, and is not affected by its compo¬ 
sition. 

2.4.4. M s heii and 4,o/i//,o - panel (f) ofFigure\3\ 

Decreasing M s h e ii decreases both the available thermonu¬ 
clear energy and the binding energy |.E b i n |. Ho wever, the ob¬ 
tained KE roughly follows |£bin| (see Section [2.4.51) . Mm 
is not sensitive to M s h e ii, unless weak explosion is obtained. 
For 4.oA//,o = 10 3 the maximal He-O shell mass is only 
Msheii.max ~ 1-02 Mq, but at a given M sbe n the results depend 
weakly on 4,o/4/,o- 

2.4.5. Mbase and phase - Figure\J\ 

In this section 4,o/f//,o = 100, Xu e = 0.5, M s h e ii = 
Msheii,max, /rot.o = 0.05, x rotj0 = 0.8, and the ranges M base £ 
[1.5,16.5] M 0 , p base £ [0.5, 2.5] • 10 4 gem -3 are scanned. 

The KE, Mm, and in-falling mass are shown in Figure [4] 
(only successful explosions are presented). The maximal pos¬ 
sible KE (Mm) increases with M base , and ranges between 
10 49 erg to 10 52 erg (negligible amounts to ss 0.5 Mq), cov¬ 
ering the observed range for the vast majority of CCSNe (see 
Section^. The in-falling mass is larger by 0.1 — 15 M 0 than 
37 ba sc ■ 

As can be seen in panel (d) of Figure [4j the KE never ex¬ 
ceeds the binding energy |£7bin| by more than a factor of 2.5. 
This can be understood by comparing the available thermonu¬ 
clear energy ~ M sbe n x MeV/ m p with the binding energy ~ 
-GMh aaff Mhp.n /rw, and noting that few x GM base /r base « 
MeV /rn p dKushnir & Katz[[2014l) . In the absence of tuning 
between the released thermonuclear energy and the min¬ 
imal KE cannot be much smaller than |£J b i n | (KE )> 0.25|E b i n | 
for all calculations). Therefore, KE ~ 4 bm for CITE, in con¬ 
trast to the neutrin o mechanism, where larger KE are obtained 
for smaller | \ dFrverll 19991 iHexer et al3l20f)l . 

3. COMPARISON TO OBSERVATIONS AND PRELIMINARY 
PREDICTIONS 

Estimates of the KE and Mn; for several observed su¬ 
pernovae are shown in Figure [3 where Type II (Type Ibc) 
supernovae, in which hydrogen is detected (not detected). 
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are marked with black crosses (x-symbols). The obse r- 


(2003^ Hendry et al. 

(2005b; 

Pastorello et al. (20051); 

Inserra etal. d20lll); 

Pastorello et al. (120121): Taddia et alJ 

(2012|); iTomasella et al. (2013 

); Dall’Ora et al. (2014b; 

iLvmanetal. (20141): Sniro et al. 

(20141): Utrobin & Chueai 


tance of 100 Mpc (to exclude rare events). The observations 
include modeling of the emitted light and are susceptible 
to systematic uncertainties. Moreover, the distribution of 
the sample in the KE-M N ; plane does not represent relative 
rates of events. Nevertheless, the range of observed KE 
(10 50 — few x 10 52 erg) and Mm (10~ 3 — 1 M 0 ), and the 
gross correlation between them (higher KE leads to higher 
Mm), can be deduced. The lower limits of the KE and Mm 
may be observationally biased. 

The calculated KE and Mm from Section [2] are compared 
to observations. Because of the sensitivity of the calculated 
resu lts to t he numerical treatment of the inner boundary (Sec- 
tion l2.3.ll) and the partial parameter scan, this comparison, as 
well as the predictions made, are preliminary. As can be seen, 
the range of observed KE and Mm, as well as the gross cor¬ 
relation between them, can be obtained from CITE. It is also 
possibl e to ob tain KE « few x 10 49 erg without tuning (see 
Section l2.4.5l) for small values of Mb ase . One feature of CITE 
is that KE scales with Mb ase , i.e. stronger explosions origi¬ 
nated from higher mass progenitors, as observations suggest 
( Ppznansfl|2013|), and in contr ast to the neutrino mechanism 


CFrverll 199 C AiHeger et alJfTOOl . 


The calculated KE and remnant mass are shown in panel 
(b) of Figure 0 The remnant mass in the neutron star (NS) 
regime, M grav , includes a correction due to the negative bind¬ 
ing energy by subtracting an assumed gravitational e nergy 
of 1.5 • lO 53 (Mg rav /M 0 ) 2 erg (11 .attimer& Yah ill 119891) . For 
gravitational mass above the maximal mass of a NS (taken 
here as 2.5 M 0 ) the remnant mass represents the baryonic 
mass of the in-falling material. Stronger explosions are pre¬ 
dicted to leave more massive remnants, in contrast to the 
neutrino mechanism (where strong explosions l eave neutron 
stars while weak explosions leave black holes; iFrverlfl999t 
iHeger et af1l2003l) . CITE predicts that NSs can only be pro¬ 
duced in weak (<; 10 51 erg) explosions (explaining, e.g., the 
low, < 10 5 °erg, KE of the Crab nebula: [Smithll2013l) . while 
strong (;> 10 51 erg) explosions must leave black holes (BHs). 
This indicates that the vast majority of observed extragalactic 
CCSNe (including SN1987A) left behind BHs. 
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Figure 1. A pre-collapse profile (density, temperature, enclosed mass, burning to free-fall time ratio, and specific angular momentum on the equatorial plane 
(jz=o)) that leads to a typical CCSN. For this profile Mb aS e — 6Mq, pbase = 10 4 gem -3 , tb,oA//,o = 100, A'ne = ATq — 0-5 and M S h e ii = AT s h e ii,max ~ 
3.60 Mq, leading to rbase ~ 4.62 • 10 9 cm, burning time at rbase of o ~ 7.9 • 10 2 s and a total mass of « 11.0 Mq (the profile below r = 2 ■ 10 8 cm has 
negligible effect on the results and is not shown). The obtained density, temperature, and enclosed mass pro files are similar to pre-collapse profiles of a 30 Mq 
star (dashed gray), calculated by Roni Waldman with the MESA stellar evolution code (Paxton et alj|201 ft) . The main difference is the existence of the He-0 
mixture. The rotation parameters are / ro t ,0 = 0-06 and :c ro t,o — 0-8 (ric^o — #rot ,0 x r base ~ 3.69 • 10 9 cm). 
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Figure 2. Dynamical evolution of the collapse for the initial conditions of Figure [2 Panels (a-c): logarithmic density maps at different times since collapse with 
isotope contours of He (red, Xn e =0.1, 0.2,0.45) and 56 Ni (black, X^i = 0.1). Panel (d): rate of thermonuclear energy production, Eburn (red), accumulated 
thermonuclear energy produced, Eb ur n (blue), and total KE of mass elements with positive radial (spherical) velocity, m (v r > 0) (black). 
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Figure 3. KE (black) and Mni (red) as a function of several numerical and physical parameters, panel (a): resolution, x-symbols are obtained from simulation 
with r re s = 10 10 cm; panel (b): ri nner ; panel (c): a: ro t,o; panel (d): / ro t, 0 , x-symbols are obtained from simulation with relaxed initial models; panel (e): 
composition; panel (f): M S h e ii, for tbo/tffQ = 10 2 (solid) andtb o/tff q = 10 3 (dashed). 
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(c) (d) 


Figure 4. Results of calculations with = 100, Xh c = Xo = 0.5, M s h e ii = M s h e ii,max, /rot,0 = 0.05, i ro t,o = 0.8, Mb ase G [1.5,16.5] Mg, and 

Phase G [0.5, 2.5] ■ 10 4 g cm' 3 (different colors represent different Phase )- Panel (a): KE; panel (b): Mn, ; panel (c): in-falling mass: panel (d): KE as a function 
of -£ bin . 


i im 




(a) 


(b) 


Figure 5. The uniform sample of Figure 0 (circles, different colors represent different Mb ase ) extended with the sensitivity calculations of Figure [3] (points). 
Panel (a): Mnj-KE. The observed sample, black crosses (x-symbols) for type II (type Ibc), is shown without error bars (given in Table [3]) for purpose of clarity. 
As can be seen, the range of observed KE and Mm, as well as the gross correlation between them, can be obtained from CITE. Panel (b): KE-remnant mass (see 
Section[3]for the definition of the remnant mass). 




















9 


Table 1 


The set of simulations examining the sensitivity of the Mb aS e 

= 6 Mq, pbase 

= 10 4 gem 

' 3 . 4,oA//,0 = 

t = 0.5, /rot.O 

= 0.05, 2rot,0 = 0-8 

calculation with A R 

= 28.6 km, 

nett = 5 ■ 

10 9 cm, r inner = 

: 5 • 10 7 cm. 

to various numerical and physical parameters. 

description 

^base 

s He-0 

-^shell 

total mass 

-^bin 

kinetic energy 

56 Ni mass 

in-falling mass 3 


[10 9 cm] 

[fc S ] 

[M e ] 

[Mq] 

[10 51 erg] 

[10 51 erg] 

o 

1 

to 

1 

[M e ] 

reference 

4.62 

7.81 

3.60 

11.0 

0.856 

1.29 

8.02 

6.50 

A R = 14.3 km 






1.40 

8.66 

6.45 

A R = 57.2 km 






1.37 

5.18 

6.61 

V res = 10 10 cm 






1.33 

8.33 

6.46 

/ inner = 100 km 






1.54 

16.2 

6.31 

dinner = 250 km 






1.40 

11.1 

6.38 

dinner = 1000 km 






1.19 

0.0147 

6.87 

^rot,0 = 0.4 






0.533 

14.3 

7.42 

*^rot,0 — 0.5 






0.545 

15.6 

7.22 

3?rot,0 — 0.6 






0.515 

13.6 

7.31 

#rot,0 — 0.7 






0.912 

10.6 

6.57 

Xrot.O = 0.9 






1.56 

3.96 

6.61 

^rot,0 = 1 






1.24 

0.329 

6.95 

^rot,0 = 1-1 






0.25 

0.0907 

9.31 

/rot. 0 = 0.015 






0.362 

0.211 

9.30 

/rot,0 — 0.02 






1.20 

6.92 

7.27 

/rot.O = 0.025 






1.46 

21.3 

6.95 

/rot.O = 0.03 






1.59 

17.8 

6.72 

/rot.O = 0.04 






1.47 

14.7 

6.55 

/rot.O = 0.06 






1.26 

5.23 

6.52 

/rot.O = 0.07 






1.16 

2.98 

6.55 

/rot,0 = 0.08 






0.972 

2.07 

6.63 

/rot.O — 0.09 






0.970 

1.46 

6.65 

o 

© 

II 

o 






0.742 

1.37 

6.86 

relaxed 






1.80 

9777 

03 

relaxed, / rot ,o = 0.1 






1.09 

1.16 

6.79 

21 He = 0.2 

TM 

TU0 

""315! 

11.4 

0.886 

!ITI1 

S7Z3 

03 

X He = 0.3 

4.61 

7.23 

3.56 

11.0 

0.839 

0.467 

8.12 

7.56 

X He = 0.4 

4.60 

7.51 

3.52 

10.9 

0.835 

1.13 

8.98 

6.61 

A' He = 0.6 

4.65 

8.13 

3.78 

11.3 

0.895 

1.11 

9.57 

6.65 

A Hc = 0.7 

4.69 

8.46 

4.06 

11.6 

0.953 

0.631 

8.74 

7.33 

A'He = 0.8 

4.75 

8.83 

4.50 

12.2 

1.04 

0.288 

6.59 

6.02 


2.71 

8.19 

0.444 

6.72 

0.173 

0.135 

3.68 

6.27 


3.23 

8.06 

0.929 

7.47 

0.295 

0.326 

5.32 

6.34 


3.58 

7.99 

1.42 

8.19 

0.404 

0.495 

8.03 

6.39 


3.86 

7.93 

1.89 

8.86 

0.506 

0.657 

10.7 

6.41 


4.09 

7.89 

2.36 

9.48 

0.602 

0.824 

9.20 

6.42 


4.30 

7.86 

2.81 

10.1 

0.694 

1.03 

9.18 

6.44 


4.48 

7.83 

3.25 

10.6 

0.783 

1.11 

8.28 

6.51 

tbfi/tfffi = 10 J 

3.01 

7.18 

0.389 

6.55 

0.157 

0.179 

2.99 

6.18 

tb,o/tff,0 = 10 3 

3.59 

7.08 

0.799 

7.11 

0.267 

0.355 

6.59 

6.24 

tb,o/tff,0 = 10 3 

3.82 

7.05 

1.02 

7.40 

0.320 

0.466 

8.66 

6.21 


a Without the binding energy correction used in panel (b) of Figure 0 
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Table 2 

The set of simulations with ii,o/t//,o = 100, A'i i e = 0.5, M s h e ii = MsheU.max, /rot,0 = 0.05, and a; rot ,o = 0.8. 


-^base 

Phase 

r base 

SHe-O 

Mihell 

total mass 

-^bin 

kinetic energy 

56 Ni mass 

in-falling mass 3 

[Mq] 

[10 4 gem -3 ] 

[10 9 cm] 

M 

\Mq] 

[Af®] 

]10 51 erg] 

[10 51 erg] 

[10- 2 M@\ 

[Mq\ 

1.5 

1.5 

1.50 

7.06 

0.540 

2.97 

0.0714 

0.0621 

- 

2.01 

1.5 

1.75 

1.36 

6.75 

0.297 

2.18 

0.0443 

0.0481 

- 

1.69 

1.5 

2 

1.25 

6.51 

0.187 

1.85 

0.0319 

0.0416 

- 

1.64 

1.5 

2.25 

1.16 

6.31 

0.130 

1.71 

0.0252 

0.0305 

- 

1.63 

1.5 

2.5 

1.09 

6.15 

0.0964 

1.64 

0.0210 

0.0207 

- 

1.63 

1.5 

2.75 

1.03 

6.01 

0.0755 

1.60 

0.0182 

0.0126 

~ 

1.63 

2.25 

1 

3.00 

7.88 

3.91 

9.06 

0.459 

0.133 

0.000347 

6.48 

2.25 

1.25 

2.63 

7.31 

2.21 

6.07 

0.282 

0.274 

0.000804 

2.72 

2.25 

1.5 

2.35 

6.91 

1.33 

4.43 

0.188 

0.305 

0.00130 

2.46 

2.25 

1.75 

2.12 

6.62 

0.842 

3.54 

0.135 

0.311 

0.00378 

2.42 

2.25 

2 

1.94 

6.39 

0.567 

3.05 

0.103 

0.255 

0.0153 

2.41 

2.25 

2.25 

1.79 

6.20 

0.403 

2.78 

0.0829 

0.211 

0.0198 

2.40 

2.25 

2.5 

1.67 

6.05 

0.300 

2.62 

0.0689 

0.167 

0.00950 

2.40 

2.25 

2.75 

1.56 

5.92 

0.232 

2.53 

0.0588 

0.118 

0.000793 

2.40 

3 

1 

30 

7.83 

4.11 

9.52 

07575 

0.444 

0.0462 

3.90 

3 

1.25 

3.12 

7.27 

2.38 

6.64 

0.370 

0.568 

0.0839 

3.27 

3 

1.5 

2.76 

6.88 

1.45 

5.10 

0.256 

0.537 

0.0981 

3.21 

3 

1.75 

2.48 

6.59 

0.931 

4.26 

0.187 

0.459 

0.229 

3.19 

3 

2 

2.25 

6.37 

0.627 

3.80 

0.144 

0.350 

0.550 

3.18 

3 

2.25 

2.06 

6.19 

0.442 

3.54 

0.115 

0.281 

0.956 

3.16 

3 

2.5 

1.90 

6.04 

0.324 

3.38 

0.0945 

0.209 

0.377 

3.16 

3 

2.75 

1.76 

5.91 

0.246 

3.28 

0.0797 

0.152 

0.135 

3.16 

4.5 

05 

5.03 

8772 

7.54 

16.0 

1.26 

0.942 

1.04 

5.86 

4.5 

1 

4.24 

7.81 

3.98 

10.3 

0.746 

1.04 

1.86 

4.97 

4.5 

1.25 

3.64 

7.25 

2.22 

7.56 

0.480 

0.872 

1.87 

4.81 

4.5 

1.5 

3.16 

6.88 

1.28 

6.17 

0.327 

0.665 

3.75 

4.75 

4.5 

1.75 

2.78 

6.60 

0.773 

5.45 

0.232 

0.477 

3.85 

4.71 

4.5 

2 

2.46 

6.39 

0.487 

5.07 

0.171 

0.312 

2.95 

4.70 

4.5 

2.25 

2.2 

6.21 

0.321 

4.87 

0.130 

0.215 

1.98 

4.68 

4.5 

2.5 

1.99 

6.07 

0.222 

4.75 

0.103 

0.143 

1.33 

4.67 

6 

0.5 

6.91 

10.6 

15.9 

30.8 

2.95 

0.763 

1.86 

20.6 

6 

0.75 

5.57 

8.71 

7.24 

16.7 

1.48 

1.74 

4.56 

7.10 

6 

1 

4.62 

7.81 

3.60 

11.0 

0.856 

1.29 

8.02 

6.50 

6 

1.25 

3.88 

7.27 

1.85 

8.43 

0.528 

0.917 

9.13 

6.30 

6 

1.5 

3.28 

6.91 

0.965 

7.20 

0.338 

0.588 

6.44 

6.26 

6 

1.75 

2.79 

6.64 

0.519 

6.61 

0.222 

0.328 

4.68 

6.22 

6 

2 

2.40 

6.43 

0.294 

6.33 

0.152 

0.205 

3.41 

6.18 

6 

2.25 

2.10 

6.27 

0.178 

6.20 

0.109 

0.113 

1.94 

6.16 

6 

2.5 

1.86 

6.12 

0.115 

6.13 

0.0817 

0.0231 

0.185 

6.18 

7.5 

0.5 

774? 

10.5 

15.7 

31.7 

3.34 

1.82 

5.75 

15.7 

7.5 

0.75 

5.94 

8.72 

6.82 

17.3 

1.65 

2.36 

13.0 

8.40 

7.5 

1 

4.83 

7.83 

3.14 

11.8 

0.915 

1.40 

14.4 

8.00 

7.5 

1.25 

3.94 

7.30 

1.43 

9.32 

0.525 

0.759 

10.6 

7.83 

7.5 

1.5 

3.20 

6.95 

0.642 

8.27 

0.303 

0.388 

7.43 

7.76 

7.5 

1.75 

2.62 

6.70 

0.298 

7.84 

0.180 

0.188 

4.77 

7.70 

7.5 

2 

2.19 

6.50 

0.151 

7.67 

0.114 

0.0813 

3.77 

7.67 

9 

0.5 

TM 

10.5 

15.4 

32.4 

3.70 

3.36 

15.1 

13.8 

9 

0.75 

6.21 

8.74 

6.33 

18.0 

1.77 

2.52 

21.2 

9.94 

9 

1 

4.92 

7.86 

2.64 

12.5 

0.928 

1.26 

17.6 

9.55 

9 

1.25 

3.85 

7.35 

1.02 

10.3 

0.478 

0.602 

9.91 

9.31 

9 

1.5 

2.98 

7.02 

0.377 

9.45 

0.241 

0.238 

7.47 

9.23 

9 

1.75 

2.35 

6.77 

0.153 

9.17 

0.130 

0.0441 

3.30 

9.19 

10.5 

0.5 

8.29 

10.6 

15.1 

33.2 

4.02 

4.24 

23.6 

14.7 

10.5 

0.75 

6.39 

8.75 

5.80 

18.6 

1.86 

2.48 

26.6 

11.6 

10.5 

1 

4.91 

7.90 

2.13 

13.3 

0.897 

1.18 

17.3 

11.1 

10.5 

1.25 

3.64 

7.41 

0.661 

11.3 

0.398 

0.451 

12.6 

10.8 

10.5 

1.5 

2.67 

7.10 

0.200 

10.7 

0.175 

0.110 

8.15 

10.7 

12 

0.5 

8.6 

10.6 

14.7 

33.9 

4.31 

4.93 

29.8 

15.5 

12 

0.75 

6.51 

8.78 

5.24 

19.3 

1.91 

2.35 

28.9 

13.2 

12 

1 

4.8 

7.94 

1.63 

14.1 

0.825 

0.880 

18.1 

12.5 

12 

1.25 

3.32 

7.49 

0.388 

12.5 

0.303 

0.225 

11.4 

12.2 

13.5 

0.5 

8.86 

10.6 

14.2 

34.7 

437 

5.48 

35.4 

16.9 

13.5 

0.75 

6.56 

8.81 

4.67 

19.9 

1.92 

2.19 

32.8 

14.6 

13.5 

1 

4.59 

8.00 

1.17 

15.0 

0.716 

0.630 

16.7 

14.0 

15 

0.5 

9.08 

10.6 

13.8 

35.4 

4.79 

5.34 

46.3 

19.0 

15 

0.75 

6.57 

8.84 

4.10 

20.6 

1.89 

2.13 

33.7 

16.0 

16.5 

0.5 

927 

10.6 

13.3 

36.1 

4.99 

5.81 

47.6 

20.3 

16.5 

0.75 

6.51 

8.87 

3.53 

21.3 

1.83 

1.71 

40.9 

17.6 


a Without the binding energy correction used in panel (b) of Figure|5] 
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Table 3 

A compilation from the literature of estimated KE and M\jj from the light curves 


Name 

Kinetic energy [10 51 erg] 

56 Ni mass [Mg] 

Type 

Reference 

Name 

Kinetic energy [10 51 erg] 

5 ®Ni mass [Mq] 

Type 

Reference 

69L 

o q+0.7 
-0.6 

0.082tS:“S 

II 

3 

73R 

9 ' 7 + 1-2 

' —0.9 

0.084ll;^ 4 

II 

3 

831 

i 

0.15 

Ibc 

3 

83N 

i 

0.15 

Ibc 

3 

84L 

i 

0.15 

Ibc 

3 

86 L 

1 q+0.5 
iO -0.3 

0.034l°;°i 8 

11 

3 

87A 

1.7 

0.075 

ii 

3 

88 A 

2 9+ 1,7 
z.Z-i 2 

0.062l°;°f 

11 

3 

89L 

i 9+O.6 

1 ‘ Z -0.5 

0 01 5+ u - uus 
u.uio _ 0 005 

ii 

3 

90E 

3.4 + 3 

0-0621°;^ 

11 

3 

91G 

i q+0.9 
±o -0.6 

0 022~*~° 008 
u.uzz _ 0 0Q6 

ii 

3 

92H 

3.1 +' 3 

0 12 Q “*" 9,053 

U.izy_o 037 

n 

3 

92ba 

1 q+0.5 

1 - 3 -0.4 

0.019±°;°°» 

ii 

3 

93J 

2.4++ 

0 1 3 + 0-02 
U-lO_o.oi 

n 

11 

941 

1 9+O.6 

1>Z -0.5 

0.08l°;°i 

Ibc 

11 

96cb 

2 1+ 1,6 
zl -0.9 

0 i2+°- 04 
u,±z -0.03 

n 

11 

97D 

0.9 

0.006 

ii 

3 

97ef 

8 

0.15 

Ibc 

3 

98A 

5.6 

0.11 

ii 

7 

98bw 

38.2+ 33 j 

0.76+1+ 

Ibc 

11 

99br 

0.6 

0 0016 - * -0 ' 0011 
u.uuio _ 0 0008 

ii 

3 

99cr 

1 0+0.8 
i - y -0.6 

0 QQ+O O 34 
u.uy_ 0> 027 

n 

3 

99dn 

7-3^;? 

0 12+ 001 
U - 1 Z _ 0.02 

Ibc 

11 

99em 

i+Sii 

0 036 + ° 009 
u.uoo _ 0 00g 

n 

1 

99em 

i 9+O.6 

1 -0.3 

0 042”*” 0 027 
u.u^Z—o 019 

ii 

3 

99ex 

3-e±?:i 

0-l®lo,04 

Ibc 

11 

99gi 

1 5 +°- 7 

n m q+o.013 
U.U1o_0 009 

ii 

3 

OOcb 

4 4+0.3 
^ -0.3 

0.083111H 

n 

1 

02 ap 

6.3l 3 2 ;8 

0 09 +0 ' 01 
u.uy_ 0 01 

Ibc 

11 

03Z 

0 245”*” 0018 
u.z^o _ 0 01g 

0 0063 + ° 0006 
u.uuoo _ 0 0006 

n 

4 

03bg 

3-s±i:i 

0 19+ 003 
u ‘ ■ Lt7 — 0.02 

ii 

11 

03gd 

1 4 + 0.3 

1 -0.3 

0 016 +001 
u.uiu _ 0 006 

n 

5 

03jd 

7 4+ 2 - 8 

1 -2.4 

0 51+ 01 
U - O1 _0.09 

Ibc 

11 

04aw 

6 6 +2 ' 3 
u,u — 3.3 

0 26~*”°' 05 
u - zo _0.04 

Ibc 

11 

04dk 

5.3^2 2 

n 9'7+0.05 

u,z ' -0.04 

Ibc 

11 

04dn 

7 1+3.5 
‘• 1 -3.6 

n 99 +O.O 4 
u -^ -0.03 

Ibc 

11 

04et 

9 q+0.3 
z,o -0.3 

0 068" 1 " 0 009 
U.U°o _ 0 009 

ii 

1 

04fe 

3.6±7| 

n q+0.05 

U ‘°— 0.05 

Ibc 

11 

04ff 

2-9±i:S 

n 99+0.04 
u -^— 0.03 

ii 

11 

04gq 

5.2+1 

0 14 + 00 7 

U-i^-o.05 

Ibc 

11 

05 az 

3.9tfr 

0.38l°;°« 

Ibc 

11 

05bf 


0 OQ +0 04 

u.uy _ 0 02 

Ibc 

11 

05cs 

0 43+ 0 03 
u, ^°— 0.03 

0.008211;°°“ 

ii 

1 

05cs 

0-16 t° 0 ° 0 l 

0 006”*~ 0 003 
u.uuo _ 0 003 

n 

2 

05hg 

2 5 + 1 - 1 

0.7611+ 

Ibc 

11 

06T 

1 9+O.6 
i ' Z -0.5 

n 1 +0.04 
u - x — 0.02 

n 

11 

06au 

3.2 

0.073 

ii 

8 

06el 

6 4+2-6 

0 16 +0 ' 03 

10 — 0.03 

n 

11 

06ep 

4 r+ 2 - 2 

0-08+02 

Ibc 

11 

06ov 

2.4 

0.127 

n 

8 

07C 

q 0+1.6 

n 9+0.05 
u,z -0.04 

Ibc 

11 

07Y 

1.91+ 

0 05 + °- 01 

Ibc 

11 

07gr 

2.9+1 

n i +0.02 
u *— 0.01 

Ibc 

11 

07od 

0.5 

0.02 

n 

6 

07ru 

1 q+6.2 

10 —7.3 

0 52”*" 0-05 

0.05 

Ibc 

11 

07uy 

10-8 til 

0 34+0 05 

u.o^_ 0 04 

Ibc 

11 

08D 

4.5t 3 ; 7 

n i +0.02 
u -0.01 

Ibc 

11 

08ax 

2.6±fi 

0 16~*” 005 
u,±o -0.04 

n 

11 

08in 

0-505l° : |t 

0 01 5+0-005 
U.U10_o 005 

ii 

1 

08in 

0 49+ 0 098 
u '^ y — 0.098 

0 012”*~ 0 005 
U.U±z_o 005 

n 

2 

09E 

0.6 

0.04 

ii 

7 

09bb 

9-211.2 

0 31+0 05 

U, °- L — 0.04 

Ibc 

11 

09bw 

0.3 

0.022 

ii 

10 

09jf 

S-Slls 

0 24 +0 03 

Ibc 

11 

1 lbm 

14+5.7 

14 -5.6 

0 7i+ 011 
u *' ± _0.09 

Ibc 

11 

lldh 

1 5+0-8 
i -°-0.7 

0 OQ + 0 01 
u.uy_ 0 01 

n 

11 

1 lhs 

1 l”*” 1 

ii -0.5 

0 04+ 0,01 

ii 

11 

12A 

0.48 

0.011 

n 

9 

12A 

0 c 52 l: >~*” 0 ' 06 

u.ozo_ 0 06 

0 016”*” 0 002 
u.uio_o 002 

ii 

1 

12 aw 

1.5 

0.06 

n 

10 

iPTF13bvn 

1 o+0.8 
i '°-0.8 

0 07“*"° 02 
uu ' - 0.02 

Ibc 

11 
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